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Abstract 

We consider N oscillators coupled by a mean held as in the Win- 
free model. The model is governed by two parameters: the coupling 
strength k and the spectrum width 7 of the frequencies of each oscil¬ 
lator centered at 1. In the uncoupled regime, k = 0, each oscillator 
possesses its own natural frequency, and the difference between the 
phases of any two oscillators grows linearly in time. In the zero-width 
regime for the spectrum, the oscillators are simultaneously in the death 
state if and only if k is above some positive value k*. We say that N 
oscillators are synchronized if the difference between any two phases 
is uniformly bounded in time. We identify a new hypothesis for the 
existence of synchronization. The domain in ( 7 , k) of synchronization 
contains {0} x [0, k*] in its closure. Moreover the domain is indepen¬ 
dent of the number of oscillators and the distribution of the frequen¬ 
cies. We show numerically that the above hypothesis is necessary for 
the existence of synchronization. 

Keywords: Coupled oscillators, Winfree Model, synchronization, desyn¬ 
chronization, periodic orbit 

1 Introduction 

In 1967 Winfree |13j proposed a model describing the synchronization of 
a population of organisms or oscillators that interact simultaneously. We 
assume that the state of each oscillator is described by a real number, the 
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phase, that corresponds to a point on the circle. We call natural frequency, 
the frequency of each oscillator, as if it were isolated from the others. The 
natural frequencies are localized inside an interval [1 — 7 ,1 + 7 ] for some 
constant 7 that we call spectrum width. The interaction of the surrounding 
oscillators on a particular oscillator has a shape independent of this oscilla¬ 
tor. The term mean field is used for this kind of interaction or coupling. We 
denote by N the number of oscillators and by k a parameter called coupling 
strength that measures the coupling between these oscillators. 

A numerical investigation and a mathematical analysis of the transitions 
between the phases show a diagram of 4 major states (see Ariaratnam and 
Strogatz m)- A phase called death state where each oscillator is at rest, 
a phase called locking state where all the oscillators possess the same non¬ 
zero frequency, a phase called incoherence state where any two oscillators 
have different frequencies, and a phase called partial death state which is a 
mixed state where a part of the oscillators are in death state and a part is in 
incoherence state. For small values of k, the Winfree model is reduced to a 
different model, the Kuramoto model that will not be studied in this article. 
For large values of k, the oscillators are in the death state. We denote by 
At* the bifurcation parameter for the death state in the case 7 = 0. 

We say that N oscillators are synchronized if maxi<jj< 7 v \ xfit) — Xj{t)\ 
is bounded from above uniformly in time. Notice that we do not consider 
the phase of the oscillators but rather their lifts. For example, in the case 
7 = 0, our definition implies that N oscillators are always synchronized, 
independently of their initial conditions, thanks to the uniqueness property 
of periodic solutions of an ODE. The existence of synchronization has been 
studied in mw for At close to 0. We identify a new hypothesis of synchro¬ 
nization that produces a domain of parameters ( 7 , At) containing {0} x [0, At*] 
in its closure. We show in addition in a numerical example that the existence 
of synchronization disappears precisely when the hypothesis is not anymore 
realized. 

The Winfree model is given by the following differential equation 

1 ^ 

Xi = UJi - K—'^P{Xj)R{Xi) ( 1 ) 

i=i 

where P and R are 27r-periodic functions, Xi{t) is the phase of the i-th 
oscillator, and X{t) = (xi(t),... ,X]\f{t)) is the global state of the system. 
Although Xi{t) represents a scalar in the torus we actually consider 

its unique continuous lift in M, that we continue to call it xfit). When N = 1 
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and cji = 1, the Winfree model reduces to the equation 

X = 1 — kP{x)R{x). (2) 

We call locking bifurcation critical parameter, the coupling critical parameter 
K* which separates the death state and the locking state in the reduced 
Winfree model: 


K* := max{K >0:1 — kP{x)R{x) >0, Vx G M}. (3) 

(Notice that k* = +oo if maxj; P{x)R{x) < 0.) We assume through out this 
work the following hypotheses: 

HI the coupling strength k is taken in the interval (0, k*), 


H2 the natural frequencies Wj are chosen in (1 — 7 ,1 + 7 ) with 7 G (0,1), 


r-27r 


H3 synchronization hypothesis: 


P{s)R'{s) 


ds > 0, Vk G (0, K*). 


/o 1 — kP{s)R{s) 

For numerical results we use a simplified version of the Winfree model 


N 

Xj = Wj - — ^ Ppixj) sin(xi), where P^(x) = 1 + cos(x + f3), (4) 

and /3 G [ 0 , 7 r]. It is important to notice that we do not make any assump¬ 
tions on the number N of oscillators; in particular we do not assume that 
N — >■ -|-oo. We do not assume either that the natural frequencies are dis¬ 
tributed according to a particular law as it is done in [T]. Our main objective 
is to exhibit a domain in the parameter space ( 7 , n) independent of N and 
any choice of the natural frequencies satisfying H2 such that the oscilla¬ 
tors stay synchronized for all time and for any initial conditions sufficiently 
packed. 

Synchronization and locking may have several meanings or definitions 
depending on the authors. We choose the following definitions. 

Definition 1 (Synchronization). We say that the oscillators Xi{t) are syn¬ 
chronized if sup]^<jj<jv |xj(f) — Xj{t)\ is bounded from above uniformly in 
time t > 0. 


Definition 2 (Periodic locking). We say that the oscillators Xj(t) are peri¬ 
odically locked to the frequency H if they are synchronized and if there exist 
27r/n-periodic functions 4'i(t) such that 


Xi{t) = Qt + 'ki(t), Vi = 1... A^, Vt > 0. 



4 


W. OUKIL, A. Kessi, Ph. Thieullen 


We notice that in the periodic locking case the average frequency admits 
the limit n as t —?• +oo and that this limit is independent of the ith oscillator. 
The fact that the limit does exist has never been addressed mathematically 
(except of course in the death state). Our main result is a partial result in 
that direction in the locking case when k is any parameter in ( 0 , At*) and 
7 G (0, r(At)) for some T : [0, At*] —A M"''. 

Theorem 3. We consider the Winfree model given by Q and satisfying 
the hypotheses H1-H3. Then, there exists an open set U in the space of 
parameters ( 7 , At) G (0,1) x (0,At*), independent of N, containing in its 
closure {0} x [0,At*], such that for every parameter ( 7 , At) G U, for every 
N >1 and every choice of natural frequencies 

1. There exists an open set invariant by the flow, and of the form, 

■= G ■ raax\xj - Xi\ < X]} 

i=i 

where : M —A (0,1) is a lir-periodic function independent of 
N. In other words, for every parameters ( 7 , At) G U, the oscillators 
are synchronized for any initial conditions Xi( 0 ) with sufficiently small 
dispersion. 

2. There exists a particular initial condition {x*{0))fLi G and a 

common frequency > 0 such that 

xUt) = vz = 1,...,IV, vt > 0, 

where ^ : M+ —a M are C^, 2tt /Q.^^,^-periodic functions, uniformly 
bounded with respect to N. In other words, the oscillators are period¬ 
ically locked for a particular initial condition. 

Notice that the simplified Winfree model Q satisfies the hypothesis H3 
when /3 = 0 as it is proved in section 

For small values of At << 1, Quinn and Strogatz [12] use the Lindstedt’s 
Method and find an approximation value of the tangent bifurcation curve 
of the synchronization sate. We mention relative results about the Winfree 
model in [HI El (H [71 |T0| . Recently Ha, Park and Ryoo studied in [3] the 
stability of the stationary state for large coupling. 
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2 Proof of theorem [3] 

Let X = (xi, ... ,xn) be a solution of 0 . We call mean //(t) and dispersion 
5{t) of X the quantities 

1 ^ 

■= ■= (5,j(t) := x,{t) - Xj{t). (5) 

We want to build a space of parameters U of the form 

U = {( 7 , k) G (0, 1) X (0, K*) : 0 < 7 < r(K)}, 

and a 27r-periodic function called dispersion curve such that, if X{t) 

is a solution of the Winfree model 0 satisfying HI - H3, 




The dispersion curve is obtained by solving a non-autonomous affine 
differential equation with periodic coefficients. The following lemma is stan¬ 
dard. 

Lemma 4. We consider the affine differential equation 


^A(s) = a- /3(s)A(s), 


( 6 ) 


where a > 0 and /3 is a , 27r-periodic function satisfying 



Then there exists a unique C'^, positive and 2 tt- periodic function solution of 
given by 



exp ( fg ^/3(u) du) — 1 


Moreover 



where fd := max(0, —j3), and f3~^ := max(0,/3). 
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Proof. The differential equation Q is a first-order non-homogeneous linear 
differential equation. The general form of solution z{t) of Q with initial 
condition z(so) £ ^ at time sq G M is given by 

rs ps ps 

z(s) = exp ( / —l3{u)du)z{so) + a / exp ( / —f3{u)du)dt. 

J Sq tJ Sq t/ t 

If z(s) is 27r-periodic, using the 27r-periodicity of /?, one obtains a unique 
z{so) given by 

Z Sn = ct —^-ns-^- 

1 -exp(//^°+ '"-I3{u)du) 

and a unique 27r-periodic solution, positive thanks to a > 0, given by 

^ exp ( -/3(u) du)dt 

1 — exp ( fQ^—/3(u} du) 

The maximum of A is obtained using the following estimates 



r2n 


= / {h — fd ){s)ds< / I3'^[s)ds< / j3'^{s)ds, 


r2TT 


r2TT 


< (/3^-/3 )is)ds+ / P {s)ds, 

Jo Jo 


r2TT 


r2n 


< / P{s)ds+ / P {s)ds. 
Jo Jo 


□ 


The following lemma shows that, if the dispersion of a solution X(t) of 
0 is a priori uniformly bounded, S{t) < D, Vt G [0, t*], then each Sij{t) is 
a sub-solution strict of an affine differential equation as in lemma 

Lemma 5. Let us assume that for some ( 7 , k) G (0,1) x (0, k*), there exist 
D > 0 and t* > 0 such that for every t G the solution X{t) of Q 

satisfies 6{t) < D. Then for every I < i, j < N, and t G [0, t*] 

^6ij < (27 -h CkD'^) - KP{iJ.)R'{n)6ij, 

^SiA < 2-f + CKD^ + CKD, 

dt ' 


where dij{t) := Xi{t) - Xj{t), /r(t) = ^ Xk{t), 


C := IIPI 


\R" 


+ Il^^llooll-R^lloo; C Il-P^lloo ||-R||oO + ||-P||oO ll-R^I 


00 1 


00 • 
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Proof. By substracting two equations Q we obtain 


dt 


N 

= Wi - Wj - — ^ P{xk) [R{xi) - Rixj)] 

^ k=l 


N 

— Wj LOj KP(^fP)R ^ ^ Ri,j,k 

^ k=l 


where = P{ht)R'{hi){xi — Xj) — P{xk)[R{xi) — R{xj)] that we bound 

from above using the following estimates 


P{xk)[R{xi) - R{xj)] = P{ii)[R{xi) - R{xj)] + Eij^k, 

^ ||-B ||oo||-^ llool^A: ^j\ ||-B ||oo||-^ ||oo-^ ; 

P{n)[R{xi) - R{xj)] = P{n)R!{ii){xi - Xj) + Eij^k, 

\Ei,j,k\ < ||P||oo||^"||oo^?^ 


and \iOi — ujj\ < 27 . 


□ 


We estimate in the following lemma the velocity of /r. We will later find 
conditions on ( 7 , k) so that ^ > 0. The constant C is dehned in lemma 

Lemma 6. Assuming the same hypotheses as in lemma\^ we have 
1 — kP{ii)R{pl) — <7 + C'k-D, VtG[0,t*]. 

Proof. We use indeed the estimates \xi — iJ.\ < D, \1 — ^ Ylk=i < 7 , 

, .. N N 

= 1- KP{fi)R{n) - 


k=l 


k=l 


1 ^ 

E - P{tk))R{x, 

i=l 


□ 


We now assume that ( 7 , k) have been chosen so that 

1 - 7 - CkD - — > 0. (7) 

Notice that the equation ^ = 1 — K'P{xi)R{xi) corresponds to 7 = 0 and 
identical initial conditions in equation 0 ^ By a small perturbation of 7 = 0, 
and a small perturbation of the initial condition a:i( 0 ) = • • • = XAr(O), one 
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obtain ^ = 1 ± e — kP{^)R{^) and ^{xi — Xj) = — kP{p)B!{ ii){xi — 

Xj). The synchronization hypothesis H3 may be interpreted as a stability 
hypothesis after changing the time variable t to s = 

We now proceed formally. Thanks to the dehnition of k* in Q , condition 
Q implies ^ > 0 on [0, t*]. We then consider the change of variable s = ^{t) 
from [0,t*] to [Sol'S*] whre sq = /.i(O) and s* = //(t*). Let r be the inverse 
function 

^ . f [so,s*] —^ [0,U] 

' \ S l-A t(s) ’ ^ 

De&ne x*{s) = Xior(s), X*{s) = (x|(s),..., x^(s)), (5*j(s) = x*(s) -x*(s). 
With respect to the new variable s, lemma[^ admits the following equivalent 
form. 


Lemma 7. We assume the same hypotheses as in lemma and the new 
condition 0. Then 

< a( 7 i 1^, D) - I3^{s)6lj{s), Vs G [sq, s*], 

where 


a(7,K,L») 


27 + CkD^ ^ (27 + CkL>^ + C/tL>)(7 + CKL>) 

X — (1 — 7 — CkD — k/k*)(1 — k/k*) ’ 

kP{s)P!{s) 

^ 1 — kP{s)R{s) 


Proof. The chain rule gives or 

fs-jil - «P(A.)Ji(M)) = fsi, + - «P(»)ii(») - Lp). 

The definition of k* implies 1 — K,P{p)R{p) > 1 — k/k* >0. From lemma 
one obtains 

{27 + CnD^) + £i;j{l - kP(.,)R{s) - 
- 1 - nPiitjRU) 

From lemma one obtains 

—fi > 1 — kP{ij)R{ijl) — 7 — CuD > 1 — 7 — CkD — k/ 

I d I ^ 27 + CkD'^ + CkD 
ds 1 —'j — CkD — k/k^ 


□ 
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Let be the unique 27r-periodic solution of the equation 


^ = a{'y,K,D) - l3^{s)A 
as 

given by lemma The following lemma gives sufficient conditions on 
( 7 , K, D) so that maxs A^^k,d('S) < D. 

Lemma 8. There exists an open set U of parameters ( 7 , k ) whose closure 
contains { 0 } x [ 0 , k*], defined in the following way 

U := {( 7 , k) G (0,1) X (0, K*) : 0 < 7 < k-D^(k)}, where 

L{k) 


D{k) := min ^1, 


2(2 + C)/{1 - k/k^) + 2C(1 + C)k/{1 - k/k^Y 


1 -exp(- 

27rK exp(fg'' fiT(s) ds) 


L{k) := 


such that, for every ( 7 , k) G U, 

1 — 7 — C'k-D(k)-> 0, and max A^ ,, £)(^)(s) < Z)(k). 

K* sG[ 0 , 27 r] 

Proof. Thanks to lemma it is enough to check that 

a{-f,K,D{K)) 


kL[k) 


< D{k). 


We have 


D{k) < 


1 — k/k* 

2k(1 + C) 


, 7 + CkD{k) < kD(k)(1 + (5)<-(1- 




1 — 7 — CkD{k) — k/k* > -11-), 




t ^‘i{‘^l + CKD{KY) CKD{K)Yi + CKD{K)) 

a( 7 , K, D{k)) < -;--h 


< kDIkY 


1 - k/k* 

2(2 TC) ^ 2C{1 + C)k 


^(1 — 

< kD{k)L{k). 


-1 — k/k* (1 — k/k*)2 . 

Definition 9. We call dispersion curve the periodic function 

^7,k(®) • ^7,k,_D(k) ('®)) ^ 


□ 
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We prove the first part of Theorem 3 The set defined in Theorem 
is obviously open and has the shape or a tubular neighborhood about the 
line (1,1,..., 1)]R with bounded convex transverse section n{X) = fiQ. We 
want to prove that is invariant by the flow of Q. We recall that U 
denotes the set of parameters ( 7 , n) defined in lemma^ 

Proof of Theorem\^ - Item 1. Let be W(0) G and 

U ■= sup{f > 0 : V 0 < < f, 6{t') < 

Assume by contradiction that < + 00 . We use the change of variable 

s = /i(f) for every s G [sqjS*]) ■so = ^(0) and s* = Then there exist 

1 < * 0 )io < X such that Since maxA..y^K < D{k) and 

1 — 7 — CkD{k) — Kj > 0, lemmaimplies 

= a(7, K, D{k)) - /3f^{s)A^^k{s*) = -^A^^k{s*). 

There exists s < s* close enough to s* such that > Aj^k{s) or 

in other words there exists t < close enough to f* such that 5iojo(^) > 
A^^i^{^{t)). We have obtained a contradiction. □ 

We now prove the second part of Theorem We show there exists 

a periodically locked solution X{t) = (xi(f),..., XAr(t)) of the system Q 

with some initial condition in that is, there exists a common rotation 

number > 0 such that Xi{t) = + 'l'i(t), Vi = I,-- - ,A where 

Tj(t) are periodic functions of period Our strategy consists in 

constructing, by fixing the mean of A( 0 ), a compact and convex transverse 
section to the closure and a continuous Poincare map '■ 

^ by waiting the first time the mean of X{t) return to 0. We 

then use Brouwer fixed point theorem to prove the existence of a fixed point 
of P'y^K- We denote by <h^(A) the flow of the equation (§. 

Lemma 10. Let he ( 7 ,^) G U where U is defined in lemma^ Define 

S.,,, = {A G : //(A) = 0} and ¥ = {!,■ ■■,!). 

Then there exist a map (the Poincare map) P^^k ■ ^'r,K ^'r,K and a 

function (the return time map) 9^^^ ■ ^-y,K 1^"'' such that, 

^0(A(x) = P,,,(A) + 27r^, VA G 

_^ mX) < __ 

1 + 7 + k||P||oo||A||oo ^ ^ 1-7-CkT»(k)-^' 
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Proof. Let be X G such that /i(X) = 0. Let be //(t) := 
and r(s) be the inverse function of /r(t) as it has been dehned in ([^ and 
Q. We prefer to write the explicit dependence on X: = /u(t) and 

Txis) = t{s). Thanks to lemma we obtain 

1 - 7 - CkD{k) - — < flx{t) < I + J + «^||-P||oo||^||oo- 

Dehne 9{X) := rx(27r). Then fLx{t)dt = 27r implies the second 

estimate of the lemma. Let be P.y^K.{X) := — 2 'it¥. Then 

n{Py^f,{X)) = nxi9{X)) -2Tt = iixo Tx{2Tt)) - 27r = 0, 

6{P^,^{X)) = 6{^^^^\X)) < A^,^iitxi9{X)) = A.,,,(27r) = A.,,,(0). 


We have shown that is a map from into itself. 


□ 


Corollary 11. The Poincare map Py,K defined in lemma 10 admits a fixed 
point X* G 


Proof. 'S'y^K, is compact and convex; P^^^ : '^'y,K ^'y,K is continuous. By 
Brouwer fixed point theorem, Py^^ admits a fixed point in X* G We 

claim that X* 0 clSy^K- Suppose by contradiction X* G (9Sy,K- Let be 


N 




S^{t) := max(i^ 




(t) and 


2 = 1 


Then there exist 1 < fo,Jo < dV, such that = Ay^K(O). As in the 

proof of item 1 of Theorem there exists fo > 0 small enough such that, 
for every 0 < f < to; joi^) < Ay^K(li*(i))- By repeating this argument 
for every 1 < zi,ji < N satisfying the equality df j_^{t) = Ay^K(M*(t)), we 
obtain for some t' > 0, d^,{t') < Ay^K(l^*(iO)' The invariance of implies 
d^{t) < Ay^K(li*(i)) for every t> t'. We have obtained a contradiction with 
the fact that 


dfi9{Xfij) = 5*(0) = Ay,K(0) = Ay,K(27r) = Ay,K(^*(6»(X*)). □ 

Proof of Theorem^ - Item 2. Corollary |11| implies the existence of a point 
X* G and a time 0* > 0 such that ‘k®*(X*) = X* + 27rlh. By the 
uniqueness property of the solutions of an ordinary differential equation, we 
have 


$®*+*(X*) = $*(X*) + 27r^, Vt > 0. 
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Let be ^(s) := = (^i(s),-- - ,^n{s)), Vs > 0. Then T is 

periodic of period 0*. We have indeed 

+ e,) = + 27r^ - 27r^4^1K = ^(s). 

L/j); 6' 5); 

Moreover the return time 0* is uniformly bounded from above with respect 
to N as in lemma ITOl □ 


3 Numerical results 


The study of organized populations has been modeled by Winfree in m in 
1967. He observed that a week coupling between independent individuals 
tends to synchronize them. An example of biological synchronization is 
given by flashes of tropical fireflies. The synchronization appears when the 
individuals emit their flashes simultaneously with the same frequency. In 
this section we analyze a simplified Winfree model given by equation Q 
with PpiyX) = 1 + cos(x + 13) and R{x) = sin(x). 

When /3 = 0, the synchronization hypothesis H3 is satisfied. Indeed 


and 




-k{PI^R + PqE!) 
1 — kPqR 


PqR! 

Jo 1 - kPqR 


r2n 


PqR 


1 — kPqR 


ds 


r2TT 


sin^fs) 


/o 1 — k( 1 + cos(s)) sin(s) 


ds > 


TT 

3' 


The numerical value of k* is ~ 0.769 and the two constants C and C of 
lemmas and are bounded from above by 3. The theoretical domain of 
parameters U for which the Winfree model is synchronized is built in lemma 
[fusing D{k) and L{k). Its size is very small, 10“^ smaller than the size of 


the numerical domain one can expect as in figure 1(a) 


In order to analyze numerically the hypothesis H3, we discuss the model 
Q for different values of (3 G [OjVr]. We use two different order parameters 
to measure the synchronization 


rxW) : = 




dx{T) := max|xi(r) - n{T)\. 

i 
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Figure[^shows the numerical domains Ufj of synchronization for three differ¬ 
ent values of /?. The domain decreases along the 7 -direction as /? increases to 
^ as we shall see. The critical parameter giving the transition to the death 
state is defined as in equation Q by k*(/ 3) := max{P^(x)i?(x) : x G 


L-i 




(a) /3 = 0 


(b) ^ = f-0.5 



(c) ^ = f - 0.25 


Figure 1: Numerical domain of synchronization of the model Q. In 
all cases we choose initial conditions randomly distributed in [—f, f] and 
N = 100 oscillators with ordered and equidistant natural frequencies uji in 
[1 — 7,1 -|- 7 ]. Up is in grey and obtained by testing the largest 7 satisfying 
\dx{T)\ < Stt. In figure [(a)| T = 1500, (3 = 0 and k*(/3) « 0.769. In figure 


(b), T = 1800, /3 = I — 0.5 and k*(/ 3) ~ 1.936. In figure (c) T = 3500, 


/3 = f - 0.25 and k*(/3) « 2.694. 


We show in figure 2(a) the variation of the critical parameter k*(/ 5) for 

0.769. Let be 


(3 G [ 0 , 7 r]; we notice that the minimum is obtained at k*(0) 

Pp{s)R'{s 


r2n 




1 — KPp{s)R{s) 


ds. 


We also show numerically in 2(a), that the gray region corresponds to the 


domain of (/3, k) satisfying (3 G [ 0 , 7 r], n G [0,k*(/ 3)], and Hi^{(3) > 0, that 
is to the domain /3 < |. As noticed by the referee, the symmetry Hk{/3) = 
— (3) implies Hk{^) = 0 for every k G [0, k*(|)]. We choose k = 0.6 
in figure 2(b) and compute the largest 7 satisfying \dx{T)\ < Ovr, that is 


the boundary of Up at k = 0.6. We notice that the numerical domain of 
synchronization is negligible if and only if (3 G [f ,7r]. 

We study in figure the variation of the order parameter rx{l3) as (3 
increases in [0, vr] at a given time T and coupling strength k. A numerical 
value rxi(3) ~ 1 suggests a very tight cluster of almost all oscillators on the 
circle, a value rx{(3) ~ 0 suggests on the contrary symmetrically distributed 
oscillators. We observe a sharp decrease of rxi(3) at /3 = | that is when 
Hk{(3) becomes negative. 
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2 (a) 


(b) 

the critical parameter of the death state transition 
k = 0.6,T = 3x1 0^ and draw the desynchroniza¬ 
tion curve, the largest 7 for which dx{T) < Sir. We use the same experiment 
parameters as in figure 


Figure 2; We plot in 
k*(/ 3). We choose in 


2 (b) 


g 


1.5 

1 

0.5 

0 


0 




tt/2 

0 

(a) 7 o = 0 


g 


1.5 

1 

0.5 

0 


0 


JWvwW*' 


■n 12 

0 


(b) 7 i = 0.011 



(c) 72 = 0.0412 


Figure 3: The plot the order parameter rxiP) at T = 3000 and n = 0.6. 
The other parameters are chosen as in figure 


Figure 1^ shows the graph of dx{t) as a function of t G [0,3 x 10^] for 
different values of 7 and /3 at a given k = 0.6. It also shows the location of 
the oscillators on the circle at the final time T = 3 x 10^. Synchronization 
has been defined by the condition dx{t) < + 00 . For /3 > |, except 

at 7 Ri 0 , t I—)■ dx{t) does not seem to be bounded: the oscillators are 
desynchronized. This observation corroborates the disappearance of the 
synchronization domain in figure when /3 > g For 7 = 0 , the oscillators 
have identical natural frequencies; the comparison principle of ODE’s in the 
periodic case implies, if maxi<jj< 7 v |xj(0) —Xj(0)| < 27r, that dx{t) < 27r for 
every t > 0: the oscillators are always synchronized. For /? < f, we choose 
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two values of 7 which are close, but above, the desynchronization curve: 
7 i > 7(7 ~ 0.25) and 72 > 7 (f — 0.5). The plots suggest a non decreasing 
function 1 1 —)■ dx{t) going to infinity with flat parts. 

4 Conclusion 

We defined synchronization as a state where the dispersion of the oscillators 
is uniformly bounded in time. For every given conpling strength k, corre¬ 
sponding to the non death state, we showed analytically the existence of 
synchronized oscillators when the hypothesis H3 is satisfied and when the 
spectrnm width 7 is sufficiently smaller than some value depending on k. 
In addition, under the same hypothesis, we showed the existence of N pe¬ 
riodic synchronized oscillators with the same strong rotation number. We 
also showed numerically, for a parametrized simplified Winfree model at a 
given coupling strength, that the hypothesis H3 seems to be necessary for 
the existence of synchronization. 
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